The First Scube3 Mutant Mouse Line with Pleiotropic Phenotypic Alterations

The vertebrate Scube (Signal peptide, CUB, and EGF-like domain-containing protein) family consists of three independent members, Scube1–3, which encode secreted cell surface-associated membrane glycoproteins. Limited information about the general function of this gene family is available, and their roles during adulthood. Here, we present the first Scube3 mutant mouse line (Scube3N294K/N294K), which clearly shows phenotypic alterations by carrying a missense mutation in exon 8, and thus contributes to our understanding of SCUBE3 functions. We performed a detailed phenotypic characterization in the German Mouse Clinic (GMC). Scube3N294K/N294K mutants showed morphological abnormalities of the skeleton, alterations of parameters relevant for bone metabolism, changes in renal function, and hearing impairments. These findings correlate with characteristics of the rare metabolic bone disorder Paget disease of bone (PDB), associated with the chromosomal region of human SCUBE3. In addition, alterations in energy metabolism, behavior, and neurological functions were detected in Scube3N294K/N294K mice. The Scube3N294K/N294K mutant mouse line may serve as a new model for further studying the effect of impaired SCUBE3 gene function.

expressed in ectodermal, endodermal, and mesodermal derivatives, as are other members of the Scube gene family (Haworth et al. 2007). Expression of these genes has been shown to be dynamic, and both reciprocal and complementary to each other (Xavier et al. 2013;Haworth et al. 2007).
Although our understanding of the function of SCUBE3 in embryonic development as well as during adulthood is still marginal, one major role appears to be in bone development and homeostasis, with another one in neurological functions. Interestingly, human SCUBE3 maps to chromosome 6p21.3, a region that has been linked to Paget disease of bone 1 (PDB1) (Fotino et al. 1977;Tilyard et al. 1982), which is characterized by focal areas of increased bone turnover (Ralston et al. 2008). SCUBE3 function is also associated with other tissues, for example, Scube3 overexpression in transgenic mice induced cardiac hypertrophy (Yang et al. 2007), and zebrafish Scube3 was recently identified as a key regulator of fast muscle development by modulating fibroblast growth factor signaling (Tu et al. 2014). Further associations of Scube3 have been reported with hedgehog signal transduction (Johnson et al. 2012), angiogenesis (Yang et al. 2013), and the immune system (Luo et al. 2012). In addition, deregulation of SCUBE3 has been found in different tumor tissues such as lung cancer (Wu et al. 2011;Zhao et al. 2013) or renal carcinomas (Morris et al. 2011).
Although SCUBE3 seems to be involved in many different organ systems and diseases, there is no suitable mouse model so far for the study of functional alterations. Recent publications on mice lacking Scube3 did not show any obvious phenotype (Xavier et al. 2010;Xavier et al. 2013). In this study, we present the first Scube3 mutant mouse line with phenotypic alterations: Scube3 N294K/N294K . The mutant mouse line carries a recessive point mutation in Scube3 and was derived from the Munich N-ethyl-N-nitrosourea (ENU) mouse mutagenesis project (MEP, Hrabé de Angelis et al. 2000;Sabrautzki et al. 2012). A systemic phenotypic characterization (Hrabé de Angelis et al. 2015) of this new mutant mouse line annotates Scube3 gene function in mice to bone metabolism and morphology, renal function, and hearing, as well as neurological and behavioral functions and energy metabolism.

MATERIALS AND METHODS
Generation of Scube3 N294K/N294K mutants ENU mutagenesis and breeding were performed as described on a pure C3HeB/FeJ (C3H) background Sabrautzki et al. 2012;Aigner et al. 2011). Briefly, C3H mice were originally purchased from the Jackson Laboratory (Bar Harbor, ME) and ENU (Serva Electrophoresis, Heidelberg, Germany) was applied in three weekly intervals by intraperitoneal injections of 90 mg/kg body weight to 10-12 wk old male mice (G0). G0 mice were mated with wild-type C3H females to produce F1 offspring. F1 males not showing any obvious phenotypic alterations were mated with wild-type C3H females to obtain the G2 generation. We either choose 6-8 female G2 mice for matings with their F1 father or performed intercross matings of G2 mice to produce at least 20 mice (G3 families). Phenotyping for dysmorphological alterations was performed according to a standardized protocol . A mutation was confirmed by showing a Mendelian distribution of expected homozygous mutant mice. The Scube3 N294K/N294K mouse line was maintained on the C3H genetic background for more than 10 generations.

Chromosomal mapping
Homozygous carriers of the G3 generation were mated to C57BL/6J (B6) wild-type mice and the progeny (F1 generation) were intercrossed. DNA was prepared from tail tips of affected offspring (F2 generation). For chromosomal mapping, a microsatellite panel for polymorphic markers between C3H and B6 was used .

Whole exome sequencing
For enrichment of exonic sequences, we used the SureSelectXT Mouse All Exon 50 Mb kit (Agilent) followed by Illumina HiSeq2000 sequencing as 100 bp paired-end runs with an average 108 · coverage (. 93% of the target being covered . 20 ·). To search for the causative variants, we compared the sequences of one Scube3 N294K/N294K mouse to one mouse from the C3HeB/FeJ background strain.
Scube3 N294K/N294K mutants, as well as 15 male and 15 female wild-type littermate controls (Scube3 WT ), was analyzed in the primary phenotyping screen of the GMC (Hrabé de Angelis et al. 2015;Gailus-Durner et al. 2005;Fuchs et al. 2012Fuchs et al. , 2011. The tests within the phenotyping pipeline and the corresponding age of the animals, as well as references for the protocols, are shown in Supplemental Material, Table S1. For further detailed investigation of observed phenotypes, secondary tests were carried out. Details for the applied protocols are listed below.

Analysis of renal function in metabolic cages
A subgroup of 48 animals (12 mutant and control mice each of both sexes) at the age of 34 wk was subjected to a renal function test using metabolic cages for single mice (Tecniplast, Buguggiate, Italy) to collect 48 hr urine samples and monitor water and urine production. During the test, mice had free access to water and pulverized food. Tests were conducted as described previously, starting with a single blood sample collection followed by urine collection over 48 hr (Fuchs et al. 2011). Urine and plasma samples were analyzed for a set of 12 clinical chemistry parameters including concentrations of sodium, potassium, and chloride (Na, K, and Cl), Ca, P i , creatinine (CREA), U, UA, and GLUC, as well as TP and ALB, as described above.

Transcriptome analysis
Transcriptome analyses from kidney samples of four Scube3 N294K/N294K and four Scube3 WT male mice at the age of 29 wk were performed following total RNA extraction (RNAeasy Midi kit, QIAGEN). Illumina Mouse Ref8 v2.0 Expression BeadChips were employed as previously described (Horsch et al. 2008;Kugler et al. 2013). Illumina Genomestudio 2011.1 was used for data normalization (cubic spline) and statistical analysis for the identification of differential gene expression was performed with SAM (Significant Analysis of Microarrays, fold change . 1.6, FDR , 6%) (Saeed et al. 2006;Tusher et al. 2001). Overrepresented functional annotations were obtained through the use of QIAGEN's Ingenuity Pathway Analysis (IPA, QIAGEN Redwood City, www.qiagen.com/ingenuity). Expression data are available at the GEO database under GSE56402.

Feces analysis
Fecal samples were collected in single caged mice for consecutive 5 d. Feces were separated from other material (bedding and food) and desiccated in a drying oven at 60°until weight constancy before bomb calorimetric combustion (IKA C7000, Staufen, Germany). Calorie content of food samples was also determined by bomb calorimetry. Energy uptake was calculated by multiplying food energy content and total amount of food consumed. Egested energy was calculated from feces energy content multiplied by the amount of egested feces (see also Rozman et al. 2014).

Statistical analysis
The Shapiro-Wilk test, as well as histograms and quantile-quantile plots, were used to assess the normality of the investigated parameters. For continuous data meeting the assumption of normality, a two-way ANOVA with a Tukey HSD post hoc test or a linear model (including body weight as an additional covariate in some cases) was performed to test genotype-sex interaction effects. In tables, normally distributed data were expressed as means 6 SDs. Parameters with a skewed distribution were analyzed using the nonparametric Wilcoxon rank sum test and were presented as median as well as 25th percentile and 75th percentile in the tables. Categorical parameters were analyzed by a Fisher's exact test and are presented as absolute numbers. For all tests, a P-value , 0.05 was used as the level of significance. Since the primary phenotyping screen in the GMC was mainly designed as a highthroughput screen for new phenotypic alterations, a correction for multiple testing of the various parameters was not performed. Data were analyzed using R software (Version 3.0.2; Foundation of Statistical Computing, Vienna, Austria).
Housing and handling of mice was according to the German Animal Welfare Act. The animal experiments were approved by the Government of Upper Bavaria (112-02, 78-06, 144-10, and 126-11).

Data availability
Data of the complete phenotypic analysis is accessible via the website of the GMC (http://tools.mouseclinic.de/phenomap/jsp/annotation/ public/phenomap.jsf), and the Scube3 N294K/N294K mutant mouse line is accessible via the European Mouse Mutant Archive (EMMA, www.infrafrontier.eu) under EMMA ID EM:10872.

RESULTS
Generation of Scube3 N294K/N294K mice The Scube3 N294K/N294K mutation was generated in the large-scale Munich ENU mutagenesis program. The mutant line was identified by screening G3 animals for morphological abnormalities resulting in reduced body size ( Figure 1A), a shorter and kinky tail, abnormal digit positioning, and an abnormal posture when hung by the tail. Scube3 N294K/N294K variants were crossed to wild-type C3HeB/FeJ (C3H) mice, and none of the offspring showed the characteristic Scube3 N294K/N294K phenotypes. However, intercrossing these heterozygous mice resulted in a fraction of about 25% of offspring that showed the characteristic Scube3 N294K/N294K phenotypes. Therefore, we considered Scube3 N294K/N294K to be a recessive mutant mouse line, and maintained the mutation on a pure C3H genetic background for more than 10 generations.

A C-A transversion at nucleotide position 882 in exon 8 of Scube3
leads to an asparagine to lysine exchange at protein position 294 (N294K). This nonsynonymous sequence variation cosegregated with the Scube3 N294K/N294K phenotype in all 15 tested mutant mice, and was not confirmed in 10 wild-type littermates and other mice of inbred C3H, BALB/c, and B6 strains. Subsequently, this mutation was confirmed in more than 100 mutant Scube3 N294K/N294K mice during maintenance breeding and excluded in the equal number of wild-type littermates.
The N294K substitution is located within the Ca-binding EGF-like domain 7 of SCUBE3, which is highly conserved between mouse, human, zebrafish, and chicken (data not shown). Using the protein prediction programs PROVEAN (http://provean.jcvi.org; Choi et al. 2012) and PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2), this sequence variation was classified as deleterious or probably damaging with a score of 0.99.

Phenotypic analysis
A cohort of 15 male and 15 female Scube3 N294K/N294K mice and the respective number of Scube3 WT mice were analyzed in the GMC in the screens for behavior, neurology, nociception, vision and eye, dysmorphology, bone and cartilage, energy metabolism, hematology, clinical chemistry, steroid metabolism, immunology, allergy, cardiovascular system, lung function, molecular phenotyping, and pathology.
Scube3 N294K/N294K mice show skeletal abnormalities and changes in bone metabolism: In addition to the aforementioned smaller body size, shorter and kinked tails, and abnormal digit positioning, X-ray analysis of the skeleton showed no obvious craniofacial abnormalities but detected malformations of the thoracic and lumbar vertebrae, and shorter femora ( Figure 1A and Table S2C) in Scube3 N294K/N294K mice that were independent of the body size reduction. Rib fusions occurred with low penetrance (3% of observed animals, data not shown). In DEXA analysis, bone mineral density (BMD) and bone mineral content (BMC) were significantly decreased (Table S2B), but the results might be confounded by body weight (Figure 1, B and C). Several parameters in pQCT-measurements supported the findings of the DEXA analysis. The effects were weaker in 12-month-old than in 9-months-old mice (File S1 and Table S3).
Ca levels at the age of 17 wk were elevated in Scube3 N294K/N294K female mice but this finding was not confirmed in older animals. P i levels were globally decreased in male mice and ALP activities in mice of both sexes were increased with effects getting weaker with age. At 1 yr of age, Scube3 N294K/N294K mice showed a significant increase of the bone resorption marker CTX-1 (Figure 2). Scube3 N294K/N294K mice have defects in bone growth, but not in embryonic patterning: For further investigation of the observed bone abnormalities, we analyzed Scube3 N294K/N294K mice at newborn and various embryonic stages by skeletal staining. Multiple hyper-ossifications in the axial and appendicular skeleton were detected. In thoracic, lumbar, and sacral vertebrae of Scube3 N294K/N294K mice, ossification center and pedicles were more closed and fused compared to those of wild-type mice ( Figure 3B, asterisks indicate the ossification center). Further, in some thoracic vertebrae, ossification centers were split in mutant mice ( Figure 3B, two asterisks in one vertebrae). In the appendicular skeleton, the coracoid processes showed hyper-ossification, and proximal and middle phalanges in digit II started to be fused ( Figure 3F, arrowhead and black asterisks). Furthermore, metacarpal, metatarsal, and tarsal were enlarged due to hyper-ossification ( Figure 3F and Figure 3H, white asterisks). These hyper-ossification phenotypes were not detected at embryonic day E15.5 around the time when endochondral ossification starts, but were detected in E17.5 Scube3 N294K/N294K embryos (data not shown), suggesting that malformations of Scube3 N294K/N294K mice are due to defects in bone growth, but not because of disturbed embryonic patterning as in vertebral segmentation disorders.
Scube3 N294K/N294K mice show symptoms pointing to disturbances in renal function: U, CREA, and K levels, as well as a-amylase activities, were increased in plasma of Scube3 N294K/N294K mice of both sexes (Table S2B). The findings were consistent for most parameters over Figure 3 Skeletal phenotypes of newborn Scube3 N294K/N294K mice. Hyper-ossification of entire vertebrae in homozygous Scube3 N294K/N294K mice (B and D) compared to wild-type mice (A and C). Asterisks in (B) indicate remarkable hyper-ossified vertebrae. Note that split ossification centers were detected in T10 and T11 (two asterisks in one vertebra). Sternebrae of Scube3 N294K/N294K mice (D) were slightly thicker and longer than ones of wild-type mice (C). Asterisks in (D) indicate hyperossified sternebrae. Bones of upper and lower limbs in homozygous mutants (F and H) were also hyper-ossified compared to wild-type (E and G). Arrowhead indicate ossified coracoid process (E). Asterisks in (F) and (H) show representative hyper-ossified metacarpals, metatarsals, phalanx, and tarsals.
the whole life period (measured at 12, 24, 36, and 52 wk of age, Figure 2). In addition, plasma glucose levels, triglyceride, and P i concentrations were decreased in both sexes, whereas total protein, albumin, and cholesterol were decreased only in male mutants (Table S2B). To follow up these findings, we analyzed renal function in metabolic cages. Calculated 24 hr CREA clearance adjusted to body weight was normal in Scube3 N294K/N294K mice, suggesting that a similar amount of plasma was filtered per g of body mass by the kidneys of mutant and control mice within 24 hr (Table 1). However, water uptake and urine production in relation to body mass (Table 1) as well as food consumption (2.85 6 0.23 vs. 2.02 6 0.35 and 3.18 6 0.57 vs. 2.64 6 0.47 g/24 hr in mutant vs. control males and females, respectively) were significantly higher in Scube3 N294K/N294K mice. Urinary concentrations of electrolytes and U were comparable, while CREA, protein (total protein and albumin), and glucose concentrations were slightly lower, and uric acid level was significantly lower, in Scube3 N294K/N294K mice (data not shown). Consequently, calculated 24 hr excretion values per 25 g body mass and fractional excretion rates were significantly increased in Scube3 N294K/N294K mice for the electrolytes Na, K, Cl, as well as Ca, U, total protein, albumin, and glucose. Uric acid and P i excretion adjusted to body mass was found to be highly variable and did not significantly differ between genotypes (Table 1).
Transcriptome analysis suggests a role of Scube3 in reabsorption and/or excretion of urine components: We performed transcriptome analysis of kidney samples from Scube3 N294K/N294K mice. Statistical analysis of gene expression patterns identified 138 differentially expressed genes functionally classified by the following overrepresented terms: cellular development, movement, and proliferation, as well as cardiovascular system function, developmental disorder, and renal and urological disease (Table 2). Additionally, literature-based research revealed several genes associated with expression in proximal tubules (Havcr1, Has2, Met, Mep1b, Mme, and Slc29a8), reabsorption of electrolytes (proximal tubule epithelium: Cldn10 and Slc12a1), and maintenance of salt/water balance of blood (renal medulla: Acta2, Adamts1, Ehd3, Kcnj15, and Umod). All these genes were downregulated in kidney.
Scube3 N294K/N294K mice have hearing deficits, alterations in the inner and middle ear, and show behavioral abnormalities: Hearing sensitivity was assessed by auditory brainstem response (ABR) to different auditory stimuli. There were significant differences at clicks and all tested frequencies; thresholds were increased in Scube3 N294K/N294K mice of both sexes, with female Scube3 N294K/N294K mice more severely affected ( Figure  4A). Further analysis of this finding revealed alterations in inner and middle ear preparations. The inner and middle ears appeared smaller and the ossicles of the middle ear had decreased size and altered shape. The incus was smaller, and the body and head of the malleus were especially decreased in size ( Figure 4B). The middle ear cavity, the bulla, part of the temporal bone, had an irregular shape as well ( Figure 4B).
In the open field test, Scube3 N294K/N294K mice showed clearly decreased locomotor and exploratory activity, as indexed by a reduction in total distance traveled and rearing activity over the course of the 20 min test. Movement velocity was also decreased. In terms of anxietyrelated behavior, Scube3 N294K/N294K mice spent less time in the center, entered into the center less often, traveled less distance within this zone, and exhibited an increased latency to enter the center ( Figure 5 and Table S2B). In the modified hole-board test, maximum and angular velocity were decreased in Scube3 N294K/N294K mice as well as the latency to engage in risk assessment by the male mutant mice (Table S2B).
In a general observation series according to the modified SHIRPA protocol, we also observed reduced locomotor activity, particularly in female Scube3 N294K/N294K mice. In addition, we found a reduction in n grip strength, which was also influenced by the body mass reduction (Table S2, B and C and Figure 5).
Scube3 N294K/N294K mice have altered function in energy metabolism: Scube3 N294K/N294K mice had significantly reduced body mass, and considerably lower absolute fat and lean mass (Table S2C). They also showed a significant shift in body composition, particularly with regards to decreased lean mass when adjusted to body mass ( Figure 6). Interestingly, the relations between body mass and fat content, as well as body mass and lean mass, were significantly different between genotypes (Table S2C), indicating a systematic effect on body composition. Scube3 N294K/N294K mice had significantly lower glucose, triglyceride, NEFA, and glycerol concentrations in plasma (Table S2B). In the glucose tolerance test, basal fasting glucose levels were significantly decreased in Scube3 N294K/N294K mice. AUC (area under the curve) values were decreased, which indicates an improved n  The ambos (incus) has a different shape and the body of the hammer (malleus) is narrowed (arrow) and the head surface decreased. Lower panel: the whole inner and middle ear was reduced in size and, besides the alterations of the ossicles (arrow pointing toward the hammer, malleus), the auditory cavity (bulla) has an irregular shape (arrows). WT, wild-type.
glucose tolerance in Scube3 N294K/N294K mice (Table S2B). We applied indirect calorimetry to investigate energy expenditure, substrate use, and locomotor activity under home-cage conditions, and found rearing behavior decreased during the 21 hr test phase. No major effects on energy turnover could be detected when VO2 was adjusted to body mass (Table  S2C). In a separate test, food consumption and gastro-intestinal functions of Scube3 N294K/N294K mice were monitored over 5 d in single caged animals. We could not detect considerable effects on food intake, energy assimilation, and the efficiency of energy extraction from food.
Scube3 N294K/N294K mice show mild alterations in cardiovascular parameters, but no hints of deficits in heart performance or conduction: Echocardiography and electrocardiography revealed several mild alterations in Scube3 N294K/N294K mice (decreased interventricular septum width, decreased left ventricular posterior wall thickness, decreased left ventricular mass, and prolonged QRS interval duration mainly in male Scube3 N294K/N294K mice, as well as decreased diastolic ventricular dimension and decreased respiration rate in female Scube3 N294K/N294K mice, Table S2A). We found a significantly lower heart weight of Scube3 N294K/N294K mice as compared to sex-matched controls (also when normalized to tibia length). However, no alterations were found in heart performance or conduction, and thus there was no clear hint for a physiologically relevant constraint of cardiovascular function due to the mutation. Further in line, the observed alterations may be confounded by the reduced body weight of Scube3 N294K/N294K mice.
Further tested organ functions are normal in Scube3 N294K/N294K mice: The analysis of immunologically relevant parameters revealed no major changes. We observed tendencies toward decreased frequency of L-selectin expressing cells within T cell populations and increased frequency of CD11b-expressing cells within the NK cell compartment (Table S2A). We did not detect any changes in hematological parameters. Eye morphology and vision were as expected for the genetic background (C3HeB/FeJ). There were no major genotype effects on prepulse inhibition response or pain perception. No differences were detected for rotarod performance. The analysis of transepidermal water loss (TEWL) from the skin of the mice did not reveal genotype-specific differences, while hair structure and coat appeared normal. There were no alterations in lung function. Scube3 N294K/N294K mice at the age of 21 wk did not present additional histopathological alterations.

DISCUSSION
We identified an ENU-derived mouse line carrying a missense mutation in Scube3 (NM_001004366.1:c.882C . A, NP_001004366.1:p. Asn294Lys), and characterized the Scube3 N294K/N294K mutant mouse in detail in the GMC (Hrabé de Angelis et al. 2015;Gailus-Durner et al. 2005;Fuchs et al. 2012). As shown previously in a Scube3 reporter mouse line, targeted replacement of exons 2 and 3 by a lacZ-cassette demonstrated early expression of the gene in craniofacial, limb, and neural tube tissues (Xavier et al. 2010). However, these mice, as well as SCUBE3 loss-of-function mice, developed normally and did not show overt phenotypic alterations (Xavier et al. 2013(Xavier et al. , 2010 whereas Scube3 N294K/N294K mice exhibited alterations in several organ systems. One reason for this discrepancy might be owed to the different genetic backgrounds of the mice in these studies. However, a more likely explanation is the redundancy of Scube family members, which might compensate the loss of SCUBE3 in the null mutant (Xavier et al. 2013), and the nature of the mutation in the Scube3 mouse model described here. The amino acid exchange in Scube3 N294K/N294K mice lies within the EGF-like domain 7 that follows the consensus sequence D/N-X-D/N-E/Q-X m -D/N Ã -X n -Y/F for Ca-binding EGF-like domains (cbEGFs, Figure 7). These domains are known to mediate proteinprotein and protein-carbohydrate interactions in a Ca-dependent manner (Downing et al. 1996). They are also structurally important for several cellular processes such as extracellular matrix architecture or specification of cell fates (Downing et al. 1996). For example, missense mutation of conserved amino acids in cbEGFs of fibrillin 1 (FBN1) are causative for Marfan syndrome due to Ca 2+ -dependent misfolding of the protein (Whiteman et al. 2007). Likewise, the amino acid exchange in the cbEGF domain 7 in Scube3 N294K/N294K mice might lead to a reduced Ca 2+ -binding capacity and/or to conformational changes, which in turn could negatively interfere with both homo-and heteromeric interaction of SCUBE proteins (Wu et al. 2004;Yang et al. 2002), posttranslational processing (Wu et al. 2011), and interaction with transforming growth factor-b and hedgehog signaling (Johnson et al. 2012;Wu et al. 2011).
Scube3 is expressed in the cartilaginous primordia of the skeleton and regions of intramembranous bone formation in the developing craniofacial region (Haworth et al. 2007). Nevertheless, besides middle ear abnormalities, Scube3 N294K/N294K mice did not show any further craniofacial abnormalities, suggesting that SCUBE3 is dispensable for craniofacial development (Xavier et al. 2013) or that the respective functions may be assigned to regions upstream of cbEGF7 of SCUBE3. However, the Scube3 N294K/N294K mutation directly affects bone morphology and bone metabolism. Scube3 N294K/N294K mice had a significant reduction in body size and weight. Interestingly, GWAS found SCUBE3 SNPs among adult and pediatric height-associated loci (Gudbjartsson et al. 2008;Zhao et al. 2010), as well as in pigs, associated with body height, body length, and rump circumference (Wang et al. 2014). The influence of the Scube3 mutation on bone metabolism was expressed by increased ALP activities and CTX-1 values in plasma, as well as by decreased BMD and BMC in DEXA and pQCT analysis. These data are supported by the finding that SCUBE3 is expressed in Figure 6 (A) Fat mass plotted against body mass (green, Scube3 N294K/N294K males; yellow, Scube3 N294K/N294K females; blue, Scube3 WT males; and red, Scube3 WT females). (B) Lean mass plotted against body mass (green, Scube3 N294K/N294K males; yellow, Scube3 N294K/N294K females; blue, Scube3 WT males; and red, Scube3 WT females). early osteoblasts and long bones (Wu et al. 2004). Interestingly, the human SCUBE3 gene is located on chromosome 6p21.3, which was discussed to be associated with Paget disease of bone 1 (PDB1) (Wu et al. 2004;Fotino et al. 1977;Tilyard et al. 1982;Good et al. 2002). Although PDB is characterized by late onset and mostly focal bone abnormalities due to increased bone turnover (Ralston et al. 2008), mouse models for genes known to be mutated in PDB, such as Tnfrsf11a for PDB2 (OMIM 603499) and Sqstm1 for PDB3 (OMIM 601530), have bone abnormalities already at birth (Dougall et al. 1999;Li et al. 2000) or not until adulthood (Kapur et al. 2004;Durán et al. 2004;Daroszewska et al. 2011). Interestingly, cases of PDB were described to be associated with raised ALP activities and deafness (Tan et al. 2014). Therefore, we suggest Scube3 N294K/N294K mice as a new model for further studies on PDB1.
SCUBE3 maps to a region that was associated with progressive bilateral hearing loss of the mid and high frequencies (DFNA31, OMIM %608645) (Snoeckx et al. 2004). The results of the ABRhearing assessment and the morphological analysis of the middle ears of Scube3 N294K/N294K animals indicate conductive hearing loss. In Scube3 N294K/N294K mice, hearing loss might be a consequence of retarded development of middle ear cavity and ossicle malformations. By contrast, Xavier et al. (2013) did not observe any gross abnormalities in the (inner) ears of Scube3 2/2 mice at the age of E17.5. Since background-specific genetic modifiers play an important role in the development of inner ears (Kiernan et al. 2007), the differences in the genetic background could be the reason for the different observations in the two studies.
The kidneys play a central role in the regulation of mineral homeostasis and chronic kidney disease (CKD) is associated with bone metabolic disease. However, the regulatory pathways involved are still not fully understood (reviewed in Hu et al. 2013;Peacock 2010;Razzaque 2011;Rowe 2012). Elevated plasma U and CREA levels found in Scube3 N294K/N294K mice hint toward possible effects on renal function. We found increased electrolyte and water excretion, but no alteration of CREA clearance. These findings could be a sign of altered renal tubular function or might be a consequence of increased electrolyte consumption due to higher food uptake causing an upregulation of excretory mechanisms. The annotations of differentially expressed genes in kidney with reabsorption of electrolytes and salt/water balance indicate effects on tubular reabsorption. Acta2  (Yang et al. 1996), and Slc39a8 (Nebert et al. 2012) are associated with electrolyte reabsorption in the proximal tubules. Together these findings hint toward a possible role of Scube3 in kidney development and function, supported by the expression of Scube3 in the developing kidney (Haworth et al. 2007), or might be related to secondary effects on the regulation of kidney function.
The human SCUBE3 protein is involved in hedgehog and TGF-b1 signaling (Xavier et al. 2013;Haworth et al. 2007;Yang et al. 2007;Johnson et al. 2012) but, so far, nothing is known about SCUBE3 signaling in this network in mice. Transcriptome profiling analysis of kidney from Scube3 N294K/N294K mice identified potential down-stream targets of SCUBE3 that are annotated to TGF-b1 signaling [C1qtnf3 (Hofmann et al. 2011 . Reduced expression of all these genes gives evidence for a supportive role of SCUBE3 within the TGF-b1 signaling pathway by binding to transforming growth factor b receptor 2 (TGFBR2) (Yang et al. 2007).
Although we did not observe structural defects of muscle tissue in the pathological examination, the repeated findings of reduced locomotor and rearing activity as well as movement velocity might be indications for functional defects, which might also be associated with an abnormal muscle energy metabolism. Indeed, very recent findings from detailed analysis of zebrafish scube3 suggest that its gene function is crucial for fast-fiber myogenesis (Tu et al. 2014). Still, decreased activity and movement velocity could also be due to the described early developmental expression of Scube3 in the neural tube (Haworth et al. 2007).
SCUBE3 expression was observed in human cultured coronary smooth muscle cells and at low levels in the heart (Wu et al. 2004), and its overexpression led to cardiac hypertrophy in transgenic mice (Yang et al. 2007). Whereas young Scube3 transgenic mice appeared Figure 7 Scube3 N294K/N294K mice are homozygous for a (C) to (A) point mutation at nucleotide position 882 that leads to an asparagine to lysine exchange at protein position 294 (N294K) in exon 8 of Scube3 (A). The mutation affects the calcium-binding EGF-like domain VII (B), which might have effects on the capabilities to form homo-or heterodimers and block TGFb/Hedgehog signaling (C), which causes phenotypic alterations in bone development and homeostasis, hearing ability, renal function, energy metabolism, neurological functions, and behavior. EGF, epidermal growth factor; TGFb, transforming growth factor b.
phenotypically normal but showed abnormal repolarization ECG patterns as well as thickening of left-ventricular septum and posterior wall thickness when analyzed at 8 months of age (Yang et al. 2007), Scube3 N294K/N294K mice showed only very mild alterations in cardiac parameters of our phenotypic analysis. Possibly, they were too young for manifestation of a severe cardiac dysfunction.
In this study, we characterized the first Scube3 mutant mouse line with phenotypic alterations, suggesting a role of Scube3 in bone metabolism and morphology, hearing, and renal function. The observed morphological abnormalities of the skeleton, impaired bone metabolism, and hearing impairments correlate with the rare metabolic bone disorder PDB, associated with the chromosomal region of human SCUBE3. Further phenotypic alterations were observed in energy metabolism parameters, in behavior, and neurological functions (Figure 7). This new mouse model will help us to better understand SCUBE3 function and diseases related to Scube3 gene mutations.